Soliton binding and low-lying singlets in frustrated odd-legged S — \ spin tubes 
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Motivated by the intriguing properties of the vanadium spin tube Na2Va07, we show that an 
effective spin-chirality model similar to that of standard Heisenberg odd-legged S = \ spin tubes 
can be derived for frustrated inter-ring couplings, but with a spin-chirality coupling constant a that 
can be arbitrarily small. Using density matrix renormalization group and analytical arguments, we 
show that, while spontaneous dimerization is always present, solitons become bound into low- lying 
singlets as a is reduced. Experimental implications for strongly frustrated tubes are discussed. 
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Spin ladders, systems which consist of a finite number 
of coupled chains, have attracted considerable attention 
recently Spin-i ladders with an even number of legs 
are expected to have a spin gap, while ladders with an 
odd number of legs behave like spin chains at low en- 
ergy. Both predictions have been largely confirmed ex- 
perimentally. Here we have taken the ladders to have 
open boundary conditions in the rung direction, which 
is the most natural definition. In comparison, spin lad- 
ders with periodic boundary conditions in the rung di- 
rection (see Fig. QJ, often referred to as spin tubes, have 
received much less attention, mostly because the prospect 
for experimental realizations was remote. Nonetheless, it 
was noticed early on that spin tubes with an odd num- 
ber of legs are not expected to behave in the same way 
as their ladder counterparts. The crucial observation is 
that the ground state of a ring with an odd number of 
sites is not just two-fold degenerate, as for a rung in an 
odd-leg ladder, but is four-fold degenerate: In addition 
to the Kramers degeneracy of the spin-^ ground state, 
there is a degeneracy due to the two possible signs of the 
ground state momentum, leading to an extra degree of 
freedom on top of the total spin, often called the chirality 
by extension of the case of a triangle. As a consequence, 
a standard L-leg spin tube (Fig. da)), defined by the 
Hamiltonian 
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where er r ; is a spin-| operator on ring r and leg I, can be 
described in the strong-ring limit ( J' <C J) by an effective 
model, valid to first order in the inter-ring coupling J' jfj, 
defined by the Hamiltonian 



N-l 



H, 



KY,Sr- S r+1 (1 + a (t+t- +1 + h.c.)) . (1) 



Here S r are the usual spin-i operators which describe 
the total spin of ring r, while t t are pseudo-spin-^ op- 



erators acting on the chirality. The parameters of the 
model are an overall coupling constant K = 4- and a 
parameter a that measures the strength of the coupling 
between spin and chirality. For ordinary spin tubes, this 
coupling is always strong: a is equal to 4 for three-leg spin 
tubes and increases with the number of legs [8|, |9| . Us- 
ing bosonization arguments , Schulz predicted that the 
ground state should be spontaneously dimerized and the 
spectrum gapped in all sectors, a prediction supported by 
further numerical and analytical work 0, 0, [13 • In ad- 
dition, the excitations have been argued to be unbound 
solitons and the gap is always a significant fraction of J'. 
All of this remarkable physics still awaits an experimental 
realization. 

In this context, the recent synthesis of Na2V30*7 0], 
whose structure may be regarded as a spin-^ nine- leg 
spin tube, has opened up new perspectives. However, 
the properties reported so far [3| do not match the prop- 
erties predicted for standard odd-legged spin tubes. In 
particular, no spin gap could be detected in zero exter- 
nal field. This might not be too suprising, however: al- 
though the overall topology of Na2Vs07 is indeed that 
of a nine-leg spin tube, the actual geometry is quite dif- 
ferent from that of Fig. [TJa). Although ab-initio calcula- 
tions have not yet reached a consensus 0, [EJ , it is likely 
that the inter-ring coupling exhibits some kind of frustra- 
tion. Since the tubes in Na2Va07 only have a C3-axis, a 
frustrated model of the type of Fig. W[b) might be more 
appropriate. 

In this Letter, we use extensive Density Matrix Renor- 
malization Group (DMRG) simulations 13] supported 
by several analytical arguments to show that inter-ring 
frustration can have dramatic consequences for the prop- 
erties of odd-legged spin tubes. In particular, we show 
that it can reduce the spin gap and bind solitons in low- 
lying singlets. This picture might resolve some of the 
puzzles of Na2V307. 

Our starting point is to notice that, as long as the 
inter-ring coupling does not break the rotational symme- 
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FIG. 1: Two models for Na2V307i (a) non-frustrated, nine- 
leg spin tube; (b) frustrated, three-leg spin tube. 



try of the tube, the effective Hamiltonian is still given by 
Eq. (yi, but with a parameter a that can take on arbi- 
trarily small values if frustration is allowed. For instance, 
for the three-leg spin tube of Fig.^b), K = (J' + 2J")/3 
and a = 4| J' - J"\/(J' + 2J"), leading to a = 2 if each 
site is coupled to two neighbors (i.e., J' — 0) 0, and 
to a = if each site is coupled to all sites of neighboring 
rings (J' = J"). 

Therefore, we concentrate on the model of Eq. Q with 
K = 1 and consider all values of a > in the following. 
In a previous DMRG study of the effective model QJ, 
Kawano and Takahashi 8] reported the finite-size scaling 
of the spin gap for the triangular tube with a = 4. Using 
White's DMRG algorithm [13| as well, we extend their 
numerical analysis to the range < a < 20. For this pur- 
pose, we classify the lowest lying excitations according to 
the quantum numbers [S z , t z ] and study the sectors [0,0], 
[0,1], [1,0] and [1,1] for open chains with up to N = 200 
sites. From an analysis of the truncation-dependence of 
the gaps, we find that convergence is reached by keep- 
ing 250 states and thus perform the calculations up to 
that limit within six finite-system sweeps. The sum of 
the discarded density-matrix eigenvalues is smaller than 
10~ 5 in all cases. 

Above a ~ 1.4, a range that includes all non- frustrated 
spin tubes, we found that a gap is indeed present, in 
agreement with Kawano and Takahashi's analysis of the 
case a = 4, but interestingly enough, we find that the 
first excitation appears in all sectors. In other words, 
spin and chirality gaps are equal above a ~ 1.4. 

The situation changes dramatically upon reducing a. 
Below a ~ 1.4, the first excitation is no longer degen- 
erate but appears only in the sector [S z =0, t z =1]. The 
first excitation is thus a chirality excitation, and the spin 
gap As and the chirality gap A r are no longer equal. To 
examine this point further, we have performed a system- 
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FIG. 2: Extrapolated excitation gaps in the four [S z , r z ] 
sectors as a function of a. Inset: Power-law fits of the small 
a results. 



atic analysis for a < 1.4 in all sectors [0,0], [0,1], [1,0] 
and [1,1], including careful finite-size scaling. The first 
excitation is always in the sector [S Z ~Q, t z — 1] and is 
non-degenerate. Below a ~ 0.5, the second excitation 
is also non-degenerate and appears in the sector [0,0]. 
In this parameter range, the first excitation that has a 
nonzero spin quantum number is the third excited state. 
This excitation manifests itself in all sectors. By follow- 
ing this excitation as a is increased, one can determine 
that it becomes the second excitation at a ~ 0.5, and 
then the first excitation at a ~ 1.4. 

The gaps corresponding to these excitations are plot- 
ted in Fig. In extracting the gaps, some care had to 
be taken regarding finite-size effects. The results were 
fitted with polynomials in l/N, where N is the length 
of the tube. Good fits could be obtained with third or- 
der polynomials. The extrapolated values lie between 
those obtained with quadratic and quartic polynomials. 
The differences between these fits were used to define the 
error-bars shown in the inset of Fig. [21 

That these gaps correspond to very different excita- 
tions is confirmed by their a-dependences. As shown in 
the inset, the results for a < 0.1 can be fitted with power 
laws of the form A oc a b with exponents b = 1.54 ± 0.06, 
1.36 ± 0.07, and 1.11 ± 0.12 for the singlet excitations in 
the sectors [0,1], [0,0], and for the first spin excitation, 
respectively. These exponents are consistent with the 
simple fractions b = 3/2, 4/3 and 1. The large error bar 
and the value significantly larger than b = 1 for the spin 
excitation is very probably a finite-size effect. Since the 
spin gap follows a linear finite-size scaling for spin tubes 
with up to 200 rings when a is very small, the extrapo- 
lations underestimate the gap, resulting in an exponent 
larger than the actual one. 

The nature of the excitations can be further explored 
by examining the nearest-neighbor expectation values 
of the spin and pseudo-spin interactions (Si ■ Si+i) and 



i+i 



h.c). In Fig. |21 these quantities are shown 
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(a) [S z =0, t z =0], ground state (b) [S z =0, t?=0], excited state 
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FIG. 3: Nearest-neighbor expectation values of spin and chi- 
rality along a chain with 50 sites in the sectors [S z — 0, t z — 
0], [0, 1] and [1, 0] for a — 0.1. The ground state (a) is dimer- 
ized. The unaltered dimer patterns and the constrictions at 
the center of the chain in (b) and (c) indicate bound soli- 
ton states, whereas the two domain walls in (d) suggest the 
presence of two unbound solitons. 



for the lowest-lying states in the important [S z , t z ] sec- 
tors. The spin and chirality degrees of freedom alternate 
synchronously. As expected, the ground state, shown in 
Fig.[3fa), is uniformly dimerized in both the spin and the 
chirality channels. 

Excitations in a dimerized spin chain can be described 
in terms of solitons, which can be viewed as domain 
walls between two dimer coverings |l4i Il5j . In a chain 
with an even number of sites, solitons always appear in 
pairs, which can either be bound or unbound. The bound 
states can also be interpreted as excited bonds, i.e., as 
spin-triplet states in a background of dimers, and would 
therefore leave the dimer pattern unchanged. In an open 
chain, a bound soliton state is located with maximum 
probability at the center of the system, whereas its un- 
bound counterpart tries to maximize both the distance 
between the solitons and the distance to the ends of the 
chain. 

The nearest-neighbor expectation values of excited 
states in the sectors [S z ~ 0, t z = 0] and [0, 1] for a = 0.1 
are shown in Fig. IHIb) and (c). The unaltered dimer 
patterns and the constrictions at the center of the chain 
indicate the presence of bound soliton states. The very 
pronounced constriction in Fig. |3fb) suggests that the 
excitation in the ground state sector is close to un- 
binding, which we find to occur in the parameter range 
0.5 < a < 0.6. For larger values of a, all expectation val- 
ues have a structure similar to that shown in Fig. |3Jd), 
representing the lowest-lying spin excitation. Here one 



can clearly identify two domain walls at around | and | 
of the chain length, suggesting the presence of two un- 
bound solitons. 

We have checked that these features are independent 
of the chain length by studying systems of up to 200 
sites. From a complete analysis of the expectation val- 
ues throughout the whole range of a, we conclude that 
lowest-lying bound states only exist in the sectors [S z = 
0, t z = 0] and [0, 1], for which the unbinding can be ob- 
served at 0.5 < a < 0.6 and 1.3 < a < 1.4, respectively. 
The lowest-lying pure spin and combined spin-chirality 
excitations are always unbound and are therefore degen- 
erate in the thermodynamic limit. 

Interestingly, the same hierarchy of states for the 
bound and unbound soliton excitations as a function of a 
can be obtained analytically by the variational approach 
of Wang |lfj for the model that includes a next-nearest- 
neighbor coupling along the legs with relative strength 
/3 = 0.5 Since it was shown by Kawano and Taka- 
hashi || that Wang's model remains in the same phase 
as the additional interaction f3 is turned off for a = 4, 
we conjecture that the two models are in the same phase 
over the whole range of a. In the same spirit, we note 
that the effective Hamiltonian of Eq. <£Q) can also be seen 
as a special case of the recently investigated spin-orbital 
models [l6l firt fisl llfll ] . with great similarities in the roles 
of chirality and orbital degrees of freedom. 

Next, we show that the scaling of the chirality gap 
A T oc a 3 / 2 can be recovered by a simple mean-field de- 
coupling of the interaction terms S r ■ S r +iT r iz T r ¥ +1 into 

(S r - S r+ i)T^ r^ +1 +S r - S r+ i(T^ t^ +1 ) - (S r - S r+ i) (r* t^ +1 ) . 

Since, according to our numerical results, the ground 
state is spontaneously dimerized, we look for a dimerized 
solution by starting with alternating expectation values 
(S r ■ S r +i) = C s - (-l) r S s and (t+t~ +1 + t~t+ +1 ) = 
C T — (— \) r 8 Tl where Ss/ T is the alternation parameter 
in the corresponding channel. The mean-field Hamilto- 
nian then describes Heisenberg and XY-chains with al- 
ternating bond strengths. Now, for Heisenberg and XY- 
chains with alternating exchange J[l + (— l) r e], the sear- 
ings of the gap and of the alternation parameters as a 
function of e are well known Up to log- 

arithmic corrections, they are given by: As oc J(e) 2//3 , 
8s oc e 1 / 3 , A r oc Je and S T oc e. The mean- field de- 
coupling then leads to a J of order one and e oc a5 T for 
the spin part, and to J oc a and e oc £5 for the chiral- 
ity. Self-consistency then requires that Ss,S T oc a 1 / 2 and 
A T oc a 3 / 2 . This last scaling is in very good agreement 
with our DMRG result for the chirality gap, for which 
no logarithmic correction could be extracted within our 
numerical accuracy. As a further check, we have also 
extracted 5s and S T as a function of a. They are not 
strictly proportional, but we believe this is due to loga- 
rithmic corrections. Indeed, the ground state energy for 
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lying chirality excitations, which lie below the spin gap, 
are expected to lead to a number of interesting conse- 
quences such as an additional low-temperature peak in 
the specific heat. We hope that these conclusions will en- 
courage further experimental investigations of Na2Vs07 
and of other odd- legged spin tubes |12| . 
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FIG. 4: Specific heat per site as a function of the temperature 
T for a = 0.2 and 1.0 calculated by exact diagonalization of 
systems with up to 8 sites. 

small a reads 
ef/ = AaS 2 s (1 + B log S s ) + 

C (a5 T ) 4 ^ 3 (logDaSr)^ 1 — aSsS T + const. 

A, B, C, and D are constants [24| independent of a. The 
alternation parameters 5 minimizing the ground state en- 
ergy obey the relation 5 T oc <5 s (log Ss + const.), leading 
to a very good fit of our numerical results (not shown). 

In order to illustrate the influence of the chirality de- 
grees of freedom, we have calculated the specific heat 
density c as a function of temperature using exact di- 
agonalization of small systems. We have used weighted 
means between systems with an even and odd number 
of sites |2^| to extrapolate the curves to the bulk limit. 
Extrapolations for systems with up to 8 sites and two par- 
ticular values of a are shown in Fig. 0] The double-peak 
structure for small a is due to the well-separated low- 
lying chirality excitations. The distinct low-temperature 
peak progressively disappears as the chirality excitations 
get closer to the spin excitations and is completely ab- 
sent at a k 1. We also note that these low- lying chirality 
excitations, which are singlets, should be detectable in 
Raman spectroscopy. Both predictions are expected to 
apply to Na2V307 if the absence of any detectable zero- 
field spin gap in that compound is indeed a consequence 
of frustration. 

In conclusion, we have shown that the low energy prop- 
erties of odd-legged spin tubes change dramatically if 
the coupling between spin and chirality degrees of free- 
dom is reduced, as would be the case for frustrated spin 
tubes. While our results confirm that excitations in non- 
frustrated or weakly frustrated spin tubes are unbound 
solitons, they show that solitons are bound into singlet 
bound states if the coupling between spin and chirality 
is sufficiently reduced by frustration. As a consequence, 
spin and chirality gaps are no longer equal, and the low- 
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